Single-qubit lasing in the strong-coupling regime 
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Motivated by recent "circuit QED" experiments we study the lasing transition and spectral prop- 
erties of single-qubit lasers. In the strong coupling, low-temperature regime quantum fluctuations 
dominate over thermal noise and strongly influence the linewidth of the laser. When the qubit and 
the resonator are detuned, amplitude and phase fluctuations of the radiation field are coupled, and 
the phase diffusion model, commonly used to describe conventional lasers, fails. We predict pro- 
nounced effects near the lasing transition, with an enhanced linewidth and non-exponential decay 
of the correlation functions. We cover a wide range of parameters by using two complementary 
approaches, one based on the Liouville equation in a Fock state basis, covering arbitrarily strong 
coupling but limited to low photon numbers, the other based on the coherent-state representation, 
covering large photon numbers but restricted to weak or intermediate coupling. 
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I. INTRODUCTION 

A strong, coherent coupling between a superconduct- 
ing qubit and an electrical resonator was first realized 
experimentally by Wallraff et al. and Chiorescu et al. 
[H, 1H . Their work stimulated substantial theoretical 
and experimental activities [Iol - fl6| devoted to the study 
of further quantum electrodynamic effects in electric cir- 
cuits. In these "circuit QED" setups a superconducting 
qubit plays the role of an artificial atom, while the radia- 
tion field is replaced by the modes of an electric resonator. 
In some of these experiments a strong enhancement of the 
resonator field and lasing were observed [3, [l5| . 

In contrast to conventional lasers where many atoms 
are coupled weakly to the light field in an optical interfer- 
ometer, in the single-qubit laser a single superconducting 
qubit is coupled strongly to the microwave field of the 
resonator circuit. Furthermore, typical circuit QED se- 
tups operate at low temperatures, where thermal noise is 
weak and quantum fluctuations, arising from the qubit- 
resonator coupling, become dominant [17H19I ]. 

In this work we study the lasing transition and the 
spectral properties of single-qubit lasers focusing on the 
regime of strong qubit-rcsonator coupling at low temper- 
atures. In this regime we find qualitatively new behavior 
of the laser linewidth. In earlier work [13, HH we observed 
that the linewidth depends in a non-monotonous way 
on the coupling strength: approaching the lasing tran- 
sition from weak coupling one observes the well-known 
linewidth narrowing. However, for stronger coupling the 
linewidth of the single-qubit laser increases again, and 
the lasing state deteriorates, a behavior not observed in 
conventional lasers. 



If the qubit and the laser are detuned from resonance 
we find an even more surprising behavior: The linewidth, 
which generally increases with detuning, is strongly en- 
hanced near the lasing transition, and for strong coupling 
or low damping rate of the resonator shows a local maxi- 
mum as a function of detuning. In this regime, the corre- 
lation functions of the resonator decay non-exponentially 
in time. As we demonstrate below, this behavior is due 
to the coupling between phase and amplitude fluctua- 
tions, which are neglected in the commonly used phase 
diffusion model 
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In order to cover a wide range of parameters we analyze 
the properties of the laser using two approaches. One is 
based on a numerical analysis of the Liouville equation 
in a basis of Fock states. It allows for arbitrarily strong 
coupling but is limited by the size of the basis which 
we can handle. This puts constraints on the thermal 
photon number and the quality factor of the resonator. 
The other is based on the coherent-state representation. 
It covers large photon numbers (and thus low damping 
rate of the resonator) but fails in the strong coupling 
regime. The results of the two approaches coincide well 
in an interesting overlapping parameter regime. 

The paper is organized as follows: In Sec. UH we de- 
scribe the model and in Sec. IHII the two approaches em- 
ployed to solve the master equation. Stationary proper- 
ties, such as the average photon number are evaluated 
in Sec. HVl In the central part of this paper, Sec. |V] we 
investigate the spectral function and the linewidth near 
the lasing transition in the frame of the two approaches. 
The consequences of amplitude fluctuations and the fail- 
ure of the phase diffusion model are discussed. We draw 
conclusions in Sec. I VII 
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II. THE MODEL 

The single-qubit laser realized by Astafiev et al. [l4j 
consists of a Cooper pair box (CPB) [23| . characterized 
by the charging energy scale E^ and the Josephson cou- 
pling energy Ej , which is coupled capacitively to a super- 
conducting coplanar waveguide resonator with frequency 
wr. In the regime typically explored in the experiments, 
only two charge states of the CPB, indicated as |0) and 
1 2) and differing by one Cooper pair, are coherently and 
(near) resonantly coupled to the cavity. In this situation, 
the system can be modelled as a two-level system (qubit) 
coupled to a harmonic oscillator with Hamiltonian Q 

H = -^{E ch f z + Ejt x ) + hdR^a + hg T z (a + a f ) . (1) 

Here we set f z = |0)(0| - |2}(2|, t x = |0) (2| + |2) (0|, while 
a and are the annihilation and creation operators of 
photons in the resonator. The qubit-resonator coupling 
strength is denoted by go- In the considered single-qubit 
laser experiments, the qubit level spacing and cjr, are in 
the range of 10 GHz (i.e., microwave frequencies) while 
the strength of the coupling go reaches 100 MHz. This is 
much stronger than usual for conventional lasers but still 
small enough to allow for the rotating-wave approxima- 
tion. Switching to the qubit 's eigenbasis we can recast 
the Hamiltonian in the Jaynes-Cummings form (2~H ]. 

Hjc = 2^ tuJ Q (Tz ^R a ^ a + % ( cr + a + f-a^) • (2) 

Here we introduce new Pauli matrices, cr xyz and a± = 
a x ± ia y . Furthermore, the the qubit's level spacing 
huiQ = sJE^+Ej and the effective coupling strength 
g = sin# go with 9 = — arctan(£'j/£ , c h) have been in- 
troduced. In the following we allow the qubit and the 
resonator to be detuned by A = wn - wr. Although 
we consider the experiment of Ref. [14[ as a motivation 
for this analysis, our model also applies to many further 
circuit QED experiments in the la sing reg ime [5(. 

As discussed in detail in Refs. [2l|, [25|, the pumping 
mechanism implemented by Astafiev et al. flij depends 
on current injection in a superconducting single-charge 
transistor and actually involves three charge states. Be- 
sides the states |0) and |2), also a one-excess-electron 
charge state is incoherently coupled to the system and 
excited via an external voltage bias. For the present dis- 
cussion it is sufficient to model this driving simply by 
adding an incoherent excitation term to the qubit's dy- 
namics. 

In this case, within the usual Markov approximation, 
for weak system-environment coupling and using the sec- 
ular approximations, we can analyze the dynamics of a 
single-qubit laser in the frame of a Bloch-Rcdfield mas- 
ter equation (26l . [27j ■ The reduced qubit-resonator den- 
sity matrix p obeys the master equation in Lindblad form 
|28[ , which in the laboratory frame reads 



Here the Liouville super-operators £r and £q account 
for the resonator's and the qubit's dissipative processes, 

CrP= — (JVth + 1) (2apa) — afap — pc^a) 

+ - iVth (2a) pa - aa)p - paa)) , (4) 

and 

r* 

£Q P = -y (<7zP<Tz - p) 

+Y (2er + pcr_ - pa-a+ - (T-a+p) . (5) 

The dissipative dynamics of the resonator depends on 
the bare damping rate k and the thermal photon number 
iVth! while the qubit's dynamics is described by excita- 
tion, relaxation, and pure dephasing with rates r-j-, 
and r* , respectively. For later use, we also introduce the 
inverse of the T\ time, Y\ = T-j- + T^, the total dephas- 
ing rate Y v = T\/2 + T* , as well as the bare population 
inversion of the qubit To = (T-f — TjJ /T±. 

As mentioned in the Introduction, in single-qubit lasers 
the lasing transition and state depend on the coupling 
strength g and the qubit-resonator detuning A. In or- 
der to characterize the lasing state in these systems we, 
therefore, study the stationary properties of the radiation 
field, such as average photon number and photon num- 
ber fluctuations as a function of these parameters. To 
fully characterize the lasing state we further investigate 
the emission spectrum of the field, 

/oo 
die- M V(t) a ( )), (6) 
-oo 

and the linewidth of the laser radiation. This allows us 
to describe the transition between the incoherent and the 
coherent state of the resonant cavity and to draw a phase 
diagram of the lasing state. 

III. METHODS 

We address the lasing state and the resonator spec- 
trum by two different methods, by a direct numerical 
diagonalization of the master equation for the density 
matrix and by a Fokker-Planck equation approach. The 
two methods cover different parameter regimes and pro- 
vide complementary physical pictures. 

A. Direct integration of the Liouville equation 
(DILE) 

After projecting on the Fock states basis, we can recast 
the master equation ([3]) in vector form 



P = -t[Hjc,p] + £qp + 4p. (3) 



p = Gp. 



(7) 
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The reduced density matrix p is arranged as a vector p, 
and G is a superoperator acting in the space of the sys- 
tem operators. This form is convenient for the numerical 
evaluation of both static and spectral properties of the 
field. 

From the stationary solution of the master equation 
Gp a = wc obtain equilibrium values such as, e.g., the 
photon number in the resonator, (n) = Tr{crap s }. From 
Eq. (JT]) wc can also derive expressions for time-dependent 
correlation functions [29j|, e.g., (a^(t)a(O)). In order to do 
so, we employ the quantum regression theorem [H, [3(| 



(a^(t)a(O)) = Tr{ai e Gt ap 8 }. 



(8) 



To proceed it is convenient to diagonalizc the superop- 
erator G and express the product ap s in terms of eigen- 
vectors of G, ap s = ^2 k CkVk- Acting on the equation 
with the exponential exp(Gt) leads to 



e Gt ap s 



J2 c ke Xkt v k , 



where Afc is the eigenvalue corresponding to the eigen- 
vector Vk . Once the eigenvalues Afc and expansion coeffi- 
cients Cfc are known, we can easily obtain the correlation 
function for all times t. 

Wc solved Eqs. ([7]) and © numerically with results to 
be presented below. For this purpose we truncate the 
Hilbcrt space of the resonator to a finite number N of 
photon number states. In this case the superoperator 
G is represented by a AN 2 x 4A 2 -matrix, growing fast 
with N, which limits the method in our case to N < 30. 
This puts constraints to the values of rates k and Ti. 
As will be shown in Sec. IIV1 the average photon number 
(n) in the resonator obeys the relation (n) < Ti/(2k). 
Hence the damping rate k should not be too small, i.e., 
K ^ ri/iV. On the other hand, for the calculation of 
average steady state values, as well as for the linewidth 
exactly at resonance (A = 0), the calculations reduce 
to solving a system of linear equations. In these cases, 
we can cover a much higher number of photon number 
states, N < 200. 



B. Fokker-Planck equation (FPE) 

Alternatively, following the route outlined in Ref. [28| , 
we can derive a Fokker-Planck equation for the single 
qubit-laser in the coherent state representation. In this 
case we represent the qubit-oscillator density matrix as 
p(t) = f d 2 a p(a,t) |a)(a|, where p(a,t) is a 2 x 2 ma- 
trix in the qubit basis states. Substituting this expansion 
into Eq. we obtain a master equation for the operator 
p{a,t). Its trace, TrQ [/5(a,i)] = P(a,t), yields the prob- 
ability for the radiation field to be in a coherent state 

|a> - ' n 

In typical experimental configurations [14| , the dynam- 
ics of the qubit are much faster than that of the res- 
onator. As a result the qubit decays on short time scales 



to a quasi-steady state, which still depends on the slowly 
varying state of the radiation field a(t). We can adia- 
batically eliminate th e q ubit's degrees of freedom via a 
projective technique [28|]. The resulting FPE contains 
a noise term arising from the qubit fluctuations. Gen- 
erally, this noise term contains all orders of derivatives 
d/da and d/da*. If the coupling is weak compared to 
the dephasing rate of the qubit, g <C r^, the noise term 
can be truncated to second order derivatives, and the 
resulting FPE in the rotating frame reads [28| 



dP 
~dt 



2 \da 
d 2 



a(l - iA/T v )C 



dada' 



(i + A 2 /rj)(i 
(JVa. + Q) 



da 2 



Ft 



7no) 



P. 



(9) 
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Here C = 2g 2 T Q /{nT v ) and n = T 1 T (p /(4g 2 ). The 
resonator-qubit coupling leads to the noise terms 

dt((a + (t)a-(0)) q - (* + ) q (a-) q 
dt((a-(t)a_(0)) q -(a_) 2 q ), (10) 

K Jo 

which can be again calculated using the quantum regres- 
sion theorem [28|, [3(|. Here (--^q denotes the average 
over the qubit in the quasi-steady state. 

To proceed we introduce polar coordinates r and tp, de- 
noting the amplitude and phase of the radiation field in 
a coherent state, respectively. In the stationary limit, 
we expect that the field distribution function is inde- 
pendent of the phase. Hence we seek for the solutions 
limf^oo P(a, t) = P s (j")- In this way we obtain for the 
steady-state distribution 



P s (r)=JVe-*« 



(11) 



normalized such that 2ir J dr r P s (r) = 1. The expression 
for $(r) along with further details of the derivation is 
given in appendix [A"l 

The correlation function of the resonator (a'(t)a(O)) 
can be calculated by introducing the conditional prob- 
ability P(r, ip, t/ro, tpo, 0) for the radiation field to be in 
the coherent state (r, tp) at time t, given that it was in the 
state (ro, ipo) at t = 0. In terms of distribution functions, 
the correlation function is given by (3~H 

(a f (t)a(0)) = J drr 2 J dtp dp dr a r\ P s (r ) 
xe i(<po-<p) P(r,(p,t/r o ,tp o ,0) 

drr 2 W(r,t), (12) 

where W(r,t) obeys the differential equation 

dW/dt = LW. (13) 

The explicit form of the differential operator L can be ob- 
tained from the FPE equation [32| and it is given in ap- 
pendix El To solve Eq. (HU) we discretize the amplitude 



4 



r and diagonalize the L matrix on the resulting lattice. 
We can then expand W(r, t) = ^ fc CkXk{ r )e~ Xkt in eigen- 
functions of L, where Xk is the eigenvalue corresponding 
to the eigenfunction Xk , and the coefficients Ck are deter- 
mined by the initial condition W(r,t = 0) = 2irrP s (r). 

Each of the two methods described so far, based on dif- 
ferent representations of the density matrix, has its ad- 
vantages and limitations. The Fock-state representation 
used in the direct integration of the Liouville equation 
is exact but, for practical reasons, can only deal with 
low photon numbers. The coherent-state representation 
and FPE approach provide insight into the distinction 
between the quantum and classical descriptions and can 
be employed in the large photon number limit. But the 
condition g/T v <C 1 must be satisfied to allow truncating 
the derivatives to second order. Hence it fails when the 
coupling is too strong. 
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IV. PHASE DIAGRAM AND STATIONARY 
PROPERTIES 

In this Section we investigate the stationary proper- 
ties of singlc-qubit lasers for different values of the qubit- 
rcsonator coupling g and detuning A. In all results pre- 
sented below we assume iVth = and keep the bare qubit 
inversion fixed at r = 0.975. The study presented here 
serves a two-fold purpose, on one hand it aims at char- 
acterizing the state of the radiation field in the different 
parameter regimes explored in the experiment flij . on 
the other hand it gives us the opportunity to compare 
the FPE and the DILE methods. 



A. Lasing transition 

In order to study the nature of the radiation field we 
first analyze the photon distribution function, P s (r), with 
results shown in Fig. [TJ For weak coupling or large de- 
tuning the system is in the thermal regime, where the 
resonator behaves like a black-body cavity and P s (r) de- 
cays monotonically as a function of r (curve c in Fig. [1} . 
Increasing the coupling strength or decreasing the de- 
tuning brings the system into a transition regime, where 
P s (r) has a maximum at a non-zero value of r, but the 
field still has some probability at r = (curve b). A 
further change of parameters pushes the system into the 
lasing regime with a photon distribution having a negli- 
gible weight at r = 0. We note that the lasing threshold 
(solid line in Fig. [T] ) estimated within the semi-classical 
approximation (2l| . 



Figure 1: (Color online) Phase diagram (and in the inset dis- 
tribution functions) in the lasing (light), non-lasing (dark) 
and the transition regimes (intermediate) as a function of 
the detuning A = ujq — ojr and coupling g between qubit 
and resonator. The distribution functions are plotted with 
K = 10~ 4 ur, g = 0.002 uj-r, and detuning A = at point a, 



0.06 Wr at point c. 



Other 

-2 



A = 0.027 ujr at point b and A = 

parameters are F v = 1.2 x 10~ 2 cjr, I\ = 1.6 x 10 " uj-r 
to = 0.975 and Nth = 0. The solid line represents the transi 
tion curve obtained from the semi-classical approximation. 
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Figure 2: (Color online) Fano factor for different values of 
detuning and coupling between qubit and resonator. Points 
a, b, and c are defined as in Fig. \T\ 
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(14) 



coincides with the boundary of the thermal regime (found 
by evaluating whether the minimum of $(r) occurs at 
r = fH). 



A further hint at the nature of the radiation field 
in the different regimes is provided by the Fano factor 
F = ((n 2 ) — (n) 2 )/(n). As one may expect, deep in the 
thermal regime the photon distribution is simply a Bosc 
distribution and, since Nth = 0, the Fano factor equals 
1. As one approaches the transition regime, due to large 
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fluctuations in the photon number, the Fano factor is 
strongly enhanced, reaching a maximum in the transi- 
tion region. Deep in the lasing regime, the field is in 
a coherent state and the Fano factor equals one again, 
F ps 1. Fig. [2] shows that there is a wide region in the 
(A, (;)-plane where F is significantly larger than 1. In this 
region, even well above the semi-classical lasing thresh- 
old, the field is not yet in a coherent state. 



B. Photon number and fluctuations 

The lasing transition is evident in the average photon 
number (n). Results are shown in Fig. [3] for varying 
coupling strength g. In the non-lasing regime, the photon 
number is small. Remarkably, even for a singlc-qubit 
laser (n) increases sharply near the threshold before it 
saturates deep in the lasing regime. Within the semi- 
classical approximation, the saturation photon number 
can be estimated as 



2k 



(15) 



in Fig. [3] shows that both approaches describe the lasing 
transition well if the bare damping rate of the resonator is 
weak (upper panel), but for a higher damping rate (lower 
panel) the FPE approach is not sufficient. 
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Figure 3: (Color online) Average photon number on reso- 
nance obtained from the DILE (solid) and the FPE approach 
(dotted). Parameters are k = 10~ 4 ojr, in panel (a) and 
k = 5 x 10~ 4 we in panel (b). The dot-dashed line indicates 
the lasing threshold (LT) estimated using Eq. (fl4)) . 



Figure 4: (Color online) Relative fluctuations as function of g 
at resonance (panel a) and as function of A at g — 0.002 tJn 
(panel b). In both panels k — ojr. 

Of interest for the following discussions are also the 
relative fluctuations of the photon number Sn/(n) and 
of the field amplitude 8r/(r) shown in Fig. 0] (For the 
parameters used, DILE and FPE coincide well.) The 
relative fluctuations, although they are rather weak in 
the lasing regime, have important effects on the spectral 
properties of the singlc-qubit laser, which we will discuss 
in the next section. 

In the non-lasing regime, the correlations between the 
qubit and the resonator are negligible and the distribu- 
tion function for the steady state can be approximated 
as P s (r) = ire- br2 /b with 
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(16) 



2(1 + t )V~ « T2+A 2 , 

which is positive in this regime. The relative fluctuations 
of the photon number are thus given by 5n/(n) = \/l + b, 
which grow with decreasing g or increasing A. The rel- 
ative amplitude fluctuations arc constant in this regime, 
namely, Sr/(r) = \/4/7r — 1 ~ 0.5, as shown in Fig. 01 



In Fig. [3] we compare the results obtained from the 
DILE and FPE methods. Both coincide for weak cou- 
pling g. They differ when the coupling becomes stronger 
since the condition g/T v <C 1 needed for the FPE ap- 
proach is violated. The comparison of the two panels 



V. LASER LINEWIDTH 

This section is devoted to the study of the emission 
spectrum of the singlc-qubit laser radiation and the de- 
pendence of its linewidth on the qubit-resonator coupling 
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and detuning. Both in the FPE and the DILE approaches 
the correlation function of the resonator can be expanded 
as (aJ(t)a(0)) — J2k a k ex P(~^kt), where — A& are the 
eigenvalues of the discretized differential operator L de- 
fined in Appendix lAlor of the superoperator G introduced 
in Section llll Al f both in the rotating frame), respectively. 
Using this expansion in Eq. ([5]) we obtain 



1 



^ + Im(A fc )] 2 + [Re(A fc )] 5 



x {Re(a fe )Re(A fe ) + Im(a fc )[w + Im(A fe )]} . (17) 

This equation serves as starting point for the calculation 
of the resonator spectrum once the eigenvalues A& and 
the coefficients are known. 

Results for two different values of the detuning are 
plotted in Fig. [5j The spectrum is characterized by a 
linewidth, which we define - even in cases where the spec- 
trum is not Lorcntzian - to be the half-width at half- 
maximum (HWHM). We note that a detuning between 
qubit and resonator shifts the emission spectrum away 
from the bare frequency of the resonator. Near res- 
onance, this frequency shift grows linearly with the de- 
tuning a mmi. 

When analyzing the emission spectrum using the two 
approaches described above we arrive at the following 
main conclusions: 

(i) The linewidth depends in a non-monotonous way 
on the coupling strength, consistent with Refs. [2(| [U 

(ii) Away from resonance (but deep in the lasing 
regime), amplitude fluctuations contribute significantly 
to the linewidth. But they are neglected in the phase 
diffusion model. 

(iii) For strong coupling or weak damping of the res- 
onator the linewidth is strongly enhanced in the tran- 
sition regime. In this case the emission spectrum is no 
longer purely Lorcntzian. 

S e (co) [arb. unit] 



A=0 

A/w s =0.02 



Hh- 

6cj 



Figure 5: (Color online) Emission spectrum S e {oj) in the ro- 
tating frame for g = 0.003 ujr and k = 5 x 10~ 4 ojr. The 
frequency shift 5uj at A = 0.02 ojr is about 0.56 k. 
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Figure 6: (Color online) Average photon number and 
linewidth as functions of the coupling strength on resonance 
with k = 10~ 4 cjr in panel (a) and k = 5x 10~ 4 urn in panel 
(b). Results obtained from the DILE and the FPE methods 
are represented by the solid and the dashed lines, respectively. 



A. Dependence on the coupling strength 



One of the central results of Refs. [2fJ, |2l| was the non- 
monotonic dependence of the linewidth on the qubit- 
rcsonator coupling strength. This is demonstrated in 
Fig. [51 where the average photon number in the resonator 
and the linewidth of the emission spectrum arc plotted as 
functions of the coupling strength g at resonance, A = 0. 
We see that, while the photon number rapidly increases 
at the lasing transition and then saturates, the linewidth 
shows a non-monotonic behavior. For increasing but still 
weak coupling we observe the linewidth narrowing which 
is typical for lasers. However, in the deep lasing regime, 
the linewidth grows again with coupling strength, leading 
to a deterioration of the lasing state. This effect is very 
pronounced for low photon numbers in the resonator and 
strong coupling. In this regime, the linewidth can even 
become larger than the bare linewidth of the resonator. 

As described in Sec. HH the DILE can be used if the 
photon number is not too high, i.e., if the damping rate 
k of the resonator is not too small. On the other hand, 
the FPE approach requires that the coupling strength 
is much smaller than the dephasing rate of the qubit, 
g < T,,. Figs. © and (O show good agreement between 
both approaches for parameters where both are valid, i.e., 
for weak coupling and low photon number. 
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B. Linewidth as function of the detuning 

At resonance (deep in the lasing regime), a satisfac- 
tory picture of the linewidth properties can be derived 
using the phase diffusion model (PDM) [22|, HH . Within 
this model the effects of amplitude fluctuations on the 
spectrum (see below) are neglected, and the linewidth 
coincides with the phase diffusion rate, which is only af- 
fected by phase fluctuations. 

The PDM turns out to be no longer valid away from 
resonance when the quantum noise dominates. As shown 
in Fig. [3 when amplitude fluctuations are taken into ac- 
count, the laser linewidth increases with growing detun- 
ing. This is obtained both in the DILE and the FPE 
approaches and in agreement with findings of Ref. [25l . In 
contrast the PDM, although working well at resonance, 
predicts the wrong slope and curvature of the linewidth 
function of detuning. 
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This form illustrates that the coupling term vanishes 
on resonance. But off-resonance in the lasing regime, 
v(r) couples amplitude fluctuations to the phase dynam- 
ics and - as shown below - qualitatively influences the 
linewidth. 
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Figure 7: (Color online) Linewidth as a function of the de- 
tuning with k = 5 x 10~ 4 o;r, g — 0.003 cjr. 

The failure of the PDM is due to the coupling of phase 
and amplitude fluctuations. In order to describe their 
effects on the linewidth, we consider the FPE in polar 
coordinates (see appendix lAl for details), 

dP(r,(p,t) k d r r „ T . „ . d 

- = Tr¥A [Nth + Qir) + Mr)]r d-r 



at 



d -i 
+ 4v(r)— + F(r)\P{r,(p,t) 



dip 

(18) 

The mixed derivative (oc d 2 /drdip) couples the dynamics 
of phase and amplitude. The coefficient v(r) is given by 



v(r) = Im [e' 2lv R(r, y)] . 



(19) 



In the limit of full population inversion, tq = 1, and 
vanishing pure dephasing, T* = 0, it reduces to 



v(r) 



Ar 2 [(5 + A 2 /Tl)n + 2r 2 ] 
4k [(l + A 2 /r2)n +r 2 ] 3 



(20) 



Figure 8: (Color online) Linewidth as functions of the detun- 
ing obtained from the DILE approach in panel (a) and from 
the FPE method in panel (b). Parameters are given in Table. 
Ill 

In the transition regime, the linewidth is enhanced and, 
as shown in Fig. [51 exhibits a peak structure for strong 
coupling or low damping of the resonator. This phe- 
nomenon arises due to the coupling of amplitude and 
phase fluctuations. For an estimate we expand v(r) of 
Eq. (|2TJ)) (neglecting pure dephasing and setting r = 1) 
near the lasing transition in the small amplitude of the 
radiation field, r 2 <c n S a.t- The result is v(r) ~ a\r 2 /n sa t 
with expansion coefficient 



(21) 



In the derivation of the above equation, we use the rela- 
tion between the coupling strength g, the damping rate 
of the resonator k, and the detuning A at the transition, 
which is given by Eq. ([14"]) . The expansion coefficient a\ 
depends on the ratio between coupling and damping rate, 



(22) 



It grows with increasing g or decreasing k, which leads to 
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Fig [81(a) ■ 

k = 5 x 10~ 4 


Fig |8](V) • 
g = 0.003 


curve I: £ = 1.125 


g = 0.003 


k = 5.0 x 10" 4 


curve II: £ = 3.125 


g = 0.005 


k = 1.8 x 10" 4 


curve III: £ = 8 


g = 0.008 


k = 7.2 x 10" 5 


curve IV: £ = 18 


ff = 0.012 


k = 3.1 x 10" 5 



Table I: Parameters in Fig. [8] Here g and K are in unit of 



gM?=0.005 
gMj=0.008 
g/w*=0.012 




AxlO 2 /^ 



an increased coupling of the amplitude fluctuations and 
hence increased linewidth. 

In the two panels of Fig. [5] we plot the linewidth as 
a function of the detuning for different g and . cis ob- 
tained in the two approaches. Since the coefficient ai 
depends on the ratio £, we choose the parameters such 
that each curve in Fig. [5] (a) has a corresponding one in 
Fig. (b) which has the same value for £, as shown in 
Table HI The linewidths represented by curves with the 
same ratio of coupling to damping rate show similar be- 
havior around the lasing transition. When £ is small (e.g. 
£ = 1.125), the linewidth grows monotonically with de- 
tuning and no peculiar features show up. For larger val- 
ues of £, the linewidth is enhanced around the transition 
regime, which becomes more pronounced as £ increases. 
Moreover, for strong coupling, e.g., for <?/wr = 0.012, the 
peak of the linewidth can be even exceed the bare value, 
k/2. In this case, the linewidth at resonance kl(A = 0) 
is already comparable to ft/2. 

In the non-lasing regime for stron g d etuning, the 
linewidth can be well approximated by 



K 

2 



g 2 r 



p2 



A^ T0 < 2' 



(23) 



approaching the bare linewidth of the resonator, ft/2, in 
the infinite detuning limit. 



C. Non-exponential decay 

Generally, as described in Eq. (fTT]) . the emission spec- 
trum S (uj) is the sum of contributions from different 
eigenvalues. We find that both in the deep lasing and 
the non-lasing regimes, a single eigenfunction, the one 
corresponding to the eigenvalue with the smallest real 
part (denoted as Ai), is sufficient to describe the spectral 
properties. In this case, the correlation function of the 
radiation field decays exponentially in time and the spec- 
trum is a Lorentzian. However, in the transition regime, 
where the linewidth shows the non-monotonic behavior, 
more than one eigenfunction contributes to the spectrum. 

As a measure of the deviation of the correlation func- 
tion from the simple exponential form, we define the rel- 
ative distance between the HWHM linewidth kl and the 
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Figure 9: (Color online) Relative distance 7 between the 
HWHM and the real part of the first eigenvalue as a func- 
tion of the detuning for different values of g with k/uj-r = 
5 x 10 -4 (DILE) in panel (a) and for different values of k 
with g/cjn = 0.003 (FPE) in panel (b). 



real part of Ai, 



7 = 



|^l-Rc [Ai]| 



(24) 



Fig. [9] shows 7 as a function of detuning for different 
values of g and ft. At resonance (deep lasing regime) and 
far off-resonance (non-lasing regime) the first eigenvalue 
is sufficient to describe the linewidth, and the decay of 
the correlation function is simply exponential. However, 
when approaching the transition regime where a peak in 
the linewidth shows up, the weights of further eigenvalues 
increase, indicating that more eigenvalues must be taken 
into account. More details are presented in Appendix IB"1 



VI. CONCLUSION 

We have studied a single-qubit laser consisting of a 
qubit coupled strongly to the resonator. The lasing be- 
havior was investigated in a parameter space spanned by 
the coupling strength g and the detuning A. Increasing 
g or decreasing A pushes the system towards the lasing 
state, characterized by an increase of the average photon 
number and a decrease of its relative fluctuations. 

In single-qubit lasers at low temperature the linewidth 
of the emission spectrum is dominated by the quantum 
noise arising from the correlation between the qubit and 
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the resonator. On resonance, at A = 0, the linewidth 
can be described by the phase diffusion model (PDM) 
for the phase fluctuations of the laser field. In this case, 
the emission spectrum is Lorentzian with linewidth sim- 
ply given by the diffusion coefficient in the Fokker-Planck 
equation. Away from resonance but still deep in the las- 
ing regime, although the amplitude fluctuations are still 
weak, their contribution to the linewidth is magnified via 
the coupling to the phase fluctuations. In this case, the 
PDM is no longer sufficient. 

With increasing coupling strength g or decreasing 
damping rate of the resonator, the coupling between the 
phase and amplitude dynamics and, hence, the contribu- 
tion of amplitude fluctuations to the linewidth becomes 
stronger. In the transition regime to the lasing state the 
correlation function of the resonator does not simply de- 
cay exponentially in time, and the Fourier transform is 
no longer Lorentzian. Their behavior is no longer gov- 
erned by a single eigenvalue Ai, distinguished from the 
other ones by having a much smaller real part. Instead 
several eigenvalues, with similar real parts, contribute to 
the spectrum. 

In our studies, we employed two complementary ap- 
proaches in order to cover a wide range of parameters. 
The DILE method in the Fock state representation is 
free from approximations on the relative strengths of 
the coupling, damping rate of the resonator, and that 
of the qubit. A practical limitation arises from the fact 
that it can deal only with low photon numbers. In the 
large-photon number limit, the FPE approach can be 
employed. The constraint on it is set by the approxi- 
mation that the coupling strength should be lower than 
the damping rate of the qubit. In an overlapping param- 
eter regime, where both methods are sufficient, results 
obtained in the two approaches show good agreement. 



Appendix A: FPE in polar coordinates 



In polar coordinates, the FPE is given by 



dP(r, ip,t) k d 



Of 



4r dr 



d ■) 
+ 4v(r)—+F(r)jP(r,(p,t) 

+ K ( D ^ r )g-^ + Mr)g- )P(r,<p,t), 



(Al) 



where u(r) and v(r) denote the real and imaginary parts 
of e~ 2llp R, respectively, and 



F{r) = r— 
dr 

+ 2r 2 



N th + Q{r) + 2u(r) 

■rS-r +4 U (r), 



D 2 (r) 
Di(r) 

X(r) 



X(r) 

N th + Q(r) - 2u(r) 
4r 2 

AC , v(r) 



2T v X{r) 

MK + A 2 ) 
T 2 K + r2) ' 



(A2) 



In the steady state we expect the distribution function 
not to depend on the phase. It can be written in the form 
of Eq. jn) with 



$(r) 



In [N th + Q(r)+2u{r)} 

r 2 [l-Cb/X(r)]+2u(r) 



+2 / dr 



r [JVth + Q(r) + 2u(r)} 



(A3) 



The correlation function of the resonator (<v(t)a(0)) 
can be obtained by introducing (32[, 



W(r,t) = / dip dip a dr a r\ P s (r ) 

xe «<po-<P) p( r ,cp,t/r , <p ,0), (A4) 
which obeys the differential equation 
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dW 



LW 



£l{[N th + Q(r) + 2u(r)]rl 
4iv(r) + F(r)}w 

K(-D 2 (r)+iDi(r))w. (A5) 



We numerically solve this problem by discretizing the 
amplitude r and diagonalize the L matrix on the lat- 
tice. We expand W(r, t) in terms of cigenfunctions of 
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1.5 r 



1.0 



0.5 



(a) 



Re[A,]/K 
Re[A 2 ]/* 
Re[A 3 ]A 



L, namely. W(r,t) = J2 n CnXn(r)e Xnt with A„ being the 
eigenvalue corresponding to the eigenfunction \ n . The 
expansion coefficients are determined by the initial con- 
dition W(r,t = 0) = 2nrP s (r). 



Appendix B: Eigenvalues contributing to the 
linewidth 



AxlO 2 Iojr 




AxW z /a> R 

Figure 10: Real parts of the eigenvalues Ai (panel (a)) and 
their weights Wi (panel (b)) as function of the detuning. The 
results are obtained from the DILE method with k = 5 x 
10~ 4 wa and g = 0.008 cjr. 



For strong coupling or weak damping of the resonator 
the linewidth is enhanced around the transition regime. 
In this case, more than one eigenvalue contribute to the 
spectrum, and it has no longer a Lorentzian form. Since 
the eigenvalues with small real parts are of importance for 
the linewidth, we consider the first three, {A^} (ordered 
with growing real part). Their weights defined as 



Wi 



(Bl) 



determine how much the corresponding cigenfunctions 
contribute to the spectrum. As shown in Figs. [TO] and 
111! when approaching the transition regime, the weight 
of the first eigenvalue, which dominates in the non-lasing 
and deep lasing regimes, decreases, while the second or 
even the third eigenvalues becomes important. 
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Figure 11: Real parts of the eigenvalues Xi (panel (a)) and 
their weights Wi (panel (b)) as function of the detuning. The 
results are obtained from the FPE method with n = 7.2 x 
10~ 5 wa and g = 0.003 cjr. 



